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I . SUMMARY 

During  the  past  several  years,  Systems,  Science  and 
Software  (S*)  personnel  have  been  actively  engaged  in  a com- 
prehensive program  involving  computer  modeling  of  the  non- 
linear processes  that  characterize  underground  nuclear  ex- 
plosions, propagation  of  the  resultant  stress  wave  through 
realistic  earth  structures  and  prediction  of  the  ground 
motion  recorded  at  teleseismic  distances  from  an  explosion. 

The  objectives  of  the  subject  project  are  to  employ  these 
modeling  and  predictive  capabilities  in  a systematic  exeunina- 
tion  of  the  effects  of  variations  in  source  and  emplacement 
parameters  on  seismic  signals  from  underground  explosions, 
and  to  investigate  methods  for  utilizing  the  general  charac- 
teristics of  seismic  waveforms  to  obtain  reliable  yield 
estimates  for  explosions. 

The  technical  phases  necessary  to  accomplish  the  ob- 
jectives of  this  project  are  as  follows: 

1.  Conduct  a systematic  theoretical  examination 
of  material,  source  and  emplacement  parameters 
which  affect  yield-magnitude  relationships 
and  compare  the  theoretical  predictions  to 
actual  observations. 

2.  Deteinnine  and  express  uncertainties  of  yield 
estimates  in  terms  of  uncertainties  in  gross 
earth  structure,  near  source  material  proper- 
ties, and  local  source  and  receiver  structure. 

Major  accomplishments  during  the  fourth  three-month 
period  of  this  project  were  realized  in  several  different 
areas  of  research.  Of  particular  importance  was  the  exer- 
cise of  our  computational  capabilities  for  the  prediction  and 
matching  of  teleseismically  recorded  body  and  surface  waves 
from  the  recent  Pahute  Mesa  explosion,  KASSERI.  Our  objective 
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in  this  experiment  was  to  model  the  near  source  (explosion) 
nonlinear  processes,  propagate  the  resultant  stress  waves 
through  realistic  earth  structures  and  finally,  generate 
synthetic  seismograms  for  comparison  with  actual  recordings 
obtained  from  five  stations  of  the  Special  Data  Collection 
System  (SDCS)  located  at  teleseismic  distances  from  the 
Nevada  Test  Site  (NTS) . An  important  additional  aspect  of 
the  waveform  matching  exercise  for  KASSERI  is  the  comparison 
of  ground  motion  for  KASSERI  with  previous  predictions  for 
the  MAST  explosion  (Barker,  et  al. , 1976). 

The  theoretical  and  observed  amplitudes  of  both  the 
b (first  half  cycle)  and  d (maximum  cycle  within  first  3 
cycles)  phases  for  KASSERI  agree  to  well  within  a factor  of 
two  at  all  five  SDCS  stations.  In  general  the  fit  is  better 
for  the  b phase  than  for  the  d phase.  A similar  result  was 
found  for  MAST  as  reported  in  Barker,  et  al . (1975).  Minor 
adjustments  of  the  upper  mantle  model  could  improve  the 
agreement  at  all  the  SDCS  stations;  a step  that  should  be 
taken  before  computing  theoretical  seismograms  for  additional 
events . 

The  relative  scaling  of  the  observed  amplitudes  for 
KASSERI  and  MAST  at  the  SDCS  sites  is  quite  consistent  (in 
terms  of  relative  yields)  except  at  RKON.  It  will  be  very 
difficult  to  explain  the  RKON  anomaly  with  present  computa- 
tional techniques  since  the  epicentral  distance  variation 
between  events  is  only  0.1  degrees.  Waveform  similarities 
at  the  SDCS  sites  WH2YK,  CPSO  and  HNME  argue  against  attribut- 
ing the  apparent  anomaly  at  RKON  to  near-source  effects. 

The  effect  of  material  strength  on  teleseismic  body 
wave  recordings,  in  particular  the  b phase,  was  examined. 

The  amplitude  of  the  b phase  increases  with  decreasing 
strength,  with  the  rate  of  increase  being  more  rapid  at 
lower  levels  of  material  strength.  An  additional  effect  is 
that  the  period  of  the  b phase  increases  rapidly  with  decreasing 
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strength  due  to  a shift  in  the  corner  frequency  of  the  source 
spectrum. 

A comparison  of  computed  and  observed  free -field 
ground  motion  was  made  for  the  PILEDRIVER  explosion,  detonated 
in  granite  at  NTS.  The  purpose  of  this  experiment  was  to  as- 
certain whether  the  observed  free  field  recordings  obtained 
at  the  time  of  the  explosion  were  affected  by  waves  reflected 
from  the  free  surface  or  from  interfaces  between  geologic 
layers.  Comparison  of  near-field  radial  motions  with  different 
structures  indicates  that  the  observed  recordings  for  Pile- 
driver  are  unaffected  by  reflected  waves  and  are  appropriate 
choices  for  normalizing  the  source  calculations. 

A check  of  the  constitutive  modeling  embodied  in  the 
SKIPPER  code  was  made  using  free-field  ground  motion  measure- 
ments for  the  GASBUGGY  explosion  detonated  in  New  Mexico. 

Fairly  complete  material  properties  data  were  available  for 
this  event.  The  computed  reduced  displacement  potential  (RDP) , 
incorporating  all  available  information,  is  in  excellent 
agreement  with  the  observed  RDP  reported  by  Perret  (1968). 

This  result  increases  our  confidence  in  source  calculations 
used  in  predicting  teleseismic  waveforms  generated  by  under- 
ground explosions. 
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II.  INTRODUCTION 


As  stated  in  the  previous  section,  the  objectives  of 
the  subject  research  project  are  to  utilize  existing  ccxp.pu- 
tational  capabilities  to  examine  the  effects  of  various 
source  and  emplacement  parameters  on  seismic  signals  from 
underground  explosions,  and  to  devise  and  evaluate  methods 
for  utilizing  the  general  characteristics  of  seismic  wave- 
forms to  obtain  reliable  estimates  of  explosion  yields.  In 
order  to  realize  these  objectives,  activity  on  this  program 
during  the  fourth  three-month  period  of  this  contract  con- 
centrated in  the  following  areas: 

1.  The  prediction  and  matching  of  teleseismic 
body  waves  generated  by  the  NTS  explosion, 

KASSERI,  and  observed  at  the  SDCS  network. 
Comparison  of  the  results  with  those  for 
the  MAST  explosion. 

2.  Exercise  of  a recently  developed  elastic  wave 
propagation  code  including  both  near-  and  far- 
field  terms,  to  check  for  contamination  of 
close-in  ground  motion  recordings  by  reflected 
waves . 

3.  Verification  of  constitutive  modeling  using 
close-in  free  field  ground  motion  data  from  the 
GASBUGGY  explosion. 

The  plan  of  the  remainder  of  this  report  is  to  present 
technical  discussions  of  each  of  these  three  research  areas, 
followed  by  a section  summarizing  the  most  important  results 
obtained  to  date. 
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III.  TECHNICAL  DISCUSSION 


3.1  TELESEISMIC  GROUND  MOTION  PREDICTIONS  FOR  THE  KASSERI 
EXPLOSION 


3.1.1  Calculation  of  the  Equivalent  Elastic  Source 

Coupling  of  the  explosive  energy  into  elastic  waves  was 
computed  using  the  one-dimensional  (spherically  symmetric) 
finite  difference  code,  SKIPPER.  Description  of  the  techni- 
que and  the  constitutive  models  may  be  found  in  Cherry,  et 
al.  (1975) . 

KASSERI  was  detonated  in  ash  flow  tuff  at  Area  20, 
Pahute  Mesa.  The  working  point  was  well  below  the  water 
table.  Standard  measurements  for  density,  grain  density,  over- 
burden density,  water  content,  saturation  and  P wave  velocity 
were  available.  It  should  be  noted  that  the  P wave  velocity 
was  not  measured  but  was  estimated  from  events  in  similar 
media. 

No  material  strength  data  were  available.  Therefore, 
we  were  forced  to  estimate  the  material  strength  from  other 
information.  A study  of  the  effect  of  plausible  strength 
variations  on  the  seismic  signal  was  carried  out  as  part  of 
the  KASSERI  investigation  and  will  be  discussed. 

The  equation  of  state  table  (pressure  as  a function 
of  specific  volume  and  specific  internal  energy)  for  the 
saturated  tuff  was  generated  using  the  TAMEOS  scheme  (Riney, 
et  al. , 1973) . In  this  scheme  grain  density  tuff  is  mixed 
with  the  appropriate  amount  of  water,  assuming  pressure 
equilibrium  between  the  mixture  and  its  components.  From  the 
bulk  modulus  (K)  obtained  from  the  TAMEOS  table  at  the  over- 
burden pressure  and  P wave  velocity  (a) , the  shear  modulus  is 
calculated  from  pa^  = K + 4y/3.  The  hydrostatic  overburden 
pressure  is  obtained  from  the  given  depth  of  burial  and  over- 
burden density. 
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All  calculations  used  an  effective  stress  law  to  ac- 
count for  the  effect  of  pore  fluid  pressure  on  the  stress 
state.  The  material  was  assumed  to  fail  in  tension;  that  is, 
a tension  failure  model  was  invoked  if  any  principle  stress 
became  tensile. 

Since  no  data  were  available  on  the  shear  strength  of 
the  KASSERI  tuff  or  similar  materials  in  Pahute  Mesa,  a 
series  of  calculations  were  made  in  which  the  material  strength 
was  varied.  All  other  material  properties  required  for  the 
calculation  (bulk  density,  grain  density,  percent  water  by 
weight,  percent  air-filled  voids,  P wave  velocity,  overburden 
pressure,  bulk  modulus  of  overburden,  shear  modulus)  were 
held  fixed  at  the  given  values. 

Previous  calculations  for  Pahute  Mesa  rhyolite  (Cherry, 
et  al . , 1975;  Barker,  et  al . , 1976)  used  a laboratory  deter- 
mined failure  envelope  for  granite.  Taking  this  as  a start- 
ing point,  we  then  perturbed  the  failure  envelope  to  values 
that  might  be  appropriate  for  the  almost  certainly  weaker 
tuff.  The  parameters  controlling  the  failure  envelope  are 
summarized  in  Table  3.1  for  the  four  calculations  completed. 

The  parameters  in  the  table  are  defined  as  follows; 


where  J'  is  the  third  deviatoric  stress  invariant,  P is  the 

3 

pressure  including  the  overburden  and  e is  the  specific 
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TABLE  3.1 

FAILURE  ENVELOPES  FOR  KASSERI  CALCULATIONS 


Calculation  Y^(kbar)  P^(kbar)  (10 ' »ergs/cm) 


R-2924 


internal  energy.  The  Y is  the  maximum  stress  difference 
or  twice  the  maximum  shear  stress.  The  meaning  of  these 
parcimeters  is  indicated  schematically  in  Figure  3.1. 

From  the  table  we  see  that  only  the  parcimeters  Y 

a 

and  Y were  changed  from  calculation  to  calculation.  The 
m 

granite  failure  surface  was  used  for  147.  Calculations  148 

and  149  are  for  the  same  failure  surface  with  the  strength 

cut  by  a factor  of  2 and  a factor  of  4.  Calculation  146 

was  done  to  see  the  effect  of  changing  the  shape  of  the  failure 

surface  while  leaving  the  high  pressure  strength  (Y  + Y ) un- 

0 m 

changed . 

The  amplitude  of  the  reduced  velocity  potential, 

|'F(a))|,  is  plotted  in  Figure  3.2  for  the  four  sources  of 
Table  3.1.  Also  shown  is  the  rhyolite  source  used  for  the 
MAST  calculations  reported  by  Barker,  ^ al . (1976)  . The 

't' (w)  is  essentially  the  far-field  displacement  spectrum.  In 
fact, 


•f 

irij^  a log  a 4'  (1  Hz)  , 

Mg  = log  y (0.05  Hz)  , 


(3.2) 


as  has  been  demonstrated  theoretically  in  a number  of  past 
reports  (e.g.,  Bache,  et  al . , 1975a,  1975b). 

In  the  following  section  we  will  be  discussing  the 
teleseismic  body  wave  that  results  from  the  source  functions 
of  Figure  3.2.  The  surface  waves  will  be  discussed  in  a 
forthcoming  technical  report.  However,  we  can  draw  some  con- 
clusions from  the  character  of  the  source  spectra  alone.  In 
Table  3.2  are  summarized  the  main  characteristics  of  the 
spectra:  The  zero  frequency  limit  (4'^),  the  peak  value  ( | 4' | ) , 
the  ratio  of  peak  to  zero  freqeuncy  limit,  the  peak  frequency 
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Figure  3.1.  Assumed  relationship  between  the  material  strength 
(Y)  and  the  hydrodynamic  component  of  stress  (P) 
for  small  e/e^^.  When  F is  used  for  the  abcissa 
the  shock  loading  path  has  a slope  of  1-2ct. 
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(fp)  and  the  width  of  the  peak  at  the  eunplitude  halfway  between 
4'^  and  | *?  | . Comparing  146  and  147,  we  see  that  the  low 
strength  portion  of  the  strength  curve  significantly  narrows 
the  peak  while  having  a relatively  minor  effect  elsewhere. 
Comparing  147,  148,  149,  as  the  strength  is  decreased  the 
spectrum  gets  larger.  Also,  the  peak  becomes  narrower  and 
moves  to  lower  frequencies.  The  MAST  source  exhibits  behavior 
between  that  of  147  and  148. 

For  MAST  and  KASSERI  the  important  parameters  controll- 
ing the  teleseismic  coupling  are  a = 4.2  km/sec,  p = 169  kbar 
for  MAST  and  a = 3.1  km/sec  and  p = 90  kbar  for  KASSERI. 

These  values  together  with  the  source  functions  in  Figure  3.2 
and  Table  3.2  give  a first  estimate  of  the  relative  size  of 
the  body  and  surface  waves  for  these  two  events.  Detailed 
comparison  of  theoretical  and  actual  seismogreuns  is  made  in 
the  following  section. 

3.1.2  Computational  Procedure  for  Ground  Motion  Predictions 

In  our  report  on  MAST  (Barker,  et  al..,  1976)  we  out- 
lined the  procedure  followed  to  predict  the  short  and  long 
period  seismograms  at  selected  SDCS  stations.  For  KASSERI 
we  follow  the  same  procedure.  Using  the  best  available 
information  for  the  rock  properties  at  the  working  point, 
we  compute  an  equivalent  elastic  source  (RDP)  for  the  explo- 
sion. This  was  discussed  in  the  previous  section.  We  also 
construct  a layered  earth  model  for  the  source  vicinity  from 
availcQjle  information.  The  remainder  of  the  calculation  is 
identical  for  the  two  events. 

In  this  section  we  discuss  the  theoretical  body  waves 
for  KASSERI  and  compare  them  to  observations.  For  compari- 
son we  also  include  the  MAST  results  from  Barker,  et  al . 

(1976) . 
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The  source  region  crustal  structure  for  KASSERI  is 
tabulated  in  Table  3.3.  The  top  two  kilometers  of  this 
structure  are  different  from  that  for  MAST;  below  that  the 
two  structures  were  taken  to  be  the  same. 

Modeling  of  the  remainder  of  the  travel  path  is  as 
described  in  our  MAST  report.  For  the  structure  at  the  re- 
ceiver an  average  crustal  model  having  little  effect  on  the 
seismogrcuns  was  used  at  all  stations.  This  model,  which  is 
tabluated  in  Table  3.4,  was  chosen  for  the  lack  of  any  better 
information.  For  the  upper  mantle  we  chose  a slightly  modi- 
fied version  of  the  Helmberger  and  Wiggins  (1971)  model  HWNE. 
The  P wave-depth  profile  for  this  model,  together  with  the 
revised  version,  HWNE-3,  is  shown  in  Figure  3.3. 

The  other  factor  to  be  selected  is  the  parameter  T/Q 
which  characterizes  the  attenuation  along  the  path.  Results 
will  be  shown  for  two  values  of  T/Q  = 1.05  and  0.95. 

3.1.4  Comparison  of  Theoretical  and  Observed  Short  Period 

Recordings 

We  chose  the  source  148  (Figure  3.2)  to  be  the  most 
appropriate  for  the  KASSERI  event.  In  the  following  section 
we  will  compare  the  seismograms  for  this  source  to  those  ob- 
tained using  the  other  three  sources  with  different  failure 
envelopes,  but  for  the  comparison  of  theory  and  observations 
made  here,  148  is  the  source.  Note  that  this  is  the  source 

with  strength  (Y  + Y ) half  that  of  granite 

0 m 

The  wave  form  comparisons  are  shown  in  Figure  3.4. 

For  each  of  the  five  SDCS  stations  at  which  theoretical 
seismogreuns  were  made,  we  show  first  the  comparison  of  theory 
and  observations  for  MAST  from  Barker,  et  al.  (1976),  then 
the  comparison  for  KASSERI.  Comparing  seismograms  for  the 
two  events  at  each  station,  we  can  see  how  much  the  signals 
vary  theoretically  and  as  observed. 


13 


R-2924 


TABLE  3.3 

SOURCE  REGION  CRUSTAL  STRUCTURE  FOR  KASSERI 


Thickness 

(km) 


a 

(km/sec) 


(km/sec) 


P 

(q/cmM 


I 


R-2924 


TABLE  3.4 


RECEIVER  REGION  CRUSTAL  STRUCTURE 


Depth 

(km) 


Thickness 

(km) 


a 

(km/sec) 
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P velocity  versus  depth  for  the  upper  mantle  model 
HVJNE-3.  This  model  is  a slightly  modified  version 
of  HWNE  (Helmberger  and  VViggins,  1971)  with  the 
differences  as  indicated  in  the  figure. 
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The  difference  between  the  two  events  are  in  the 
source  material,  as  mentioned,  and  in  the  yield  which  is 
about  twice  as  large  for  KASSERI  as  MAST.  The  distance  to 
the  stations  of  Figure  3.4  is  approximately  10  km  greater 
for  KASSERI  than  MAST,  a seemingly  trivial  change.  The  depths 
are  H = 0.91  km  for  MAST  and  H = 1.26  km  for  KASSERI.  The 
theoretical  P-pP  lag  time  is  then  approximately  0.50  seconds 
for  MAST  and  0.82  seconds  for  KASSERI. 

Excimining  the  observations,  we  see  an  extraordinary 
change  in  the  shape  of  the  RKON  recording  between  the  two 
events.  At  the  other  stations  the  difference  are  much 
less  pronounced.  In  fact,  at  the  other  stations  one  could 
argue  that  the  waveform  differences  are  primarily  the  result 
of  the  change  in  P-pP  lag  time. 

The  theoretical  seismogrcuns  were  computed  with  T/Q  = 
0.95.  We  also  computed  the  same  stations  with  T/Q  = 1.05 
and  will  be  discussing  measurements  made  from  both  sets.  For 
a visual  comparison  as  in  Figure  3.4,  the  difference  between 
the  two  is  hardly  apparent.  At  each  station  the  instrument 
response  was  taken  from  calibration  data  provided  by  the  Pro- 
ject Officer.  The  calibration  was  done  just  prior  to  MAST 
but  was  used  for  KASSERI  as  well. 

Regarding  the  comparison  of  theoretical  and  observed 
waveforms,  we  can  only  make  the  same  comments  as  in  the  MAST 
report.  The  phenomena  that  seems  to  dominate  the  response 
at  the  nearer  four  stations  is  not  accurately  represented 
in  our  models.  Our  theoretical  signals  also  fail  to  include 
much  of  the  high  frequency  energy  that  appears  after  the 
first  few  seconds  of  the  record. 

The  important  fact  is  that  it  does  appear  that  we  have 
properly  included  the  factors  controlling  the  amplitude  of 
the  b phase,  the  first  cycle  on  the  record.  This  phase  is 
least  affected  by  phenomena  not  related  to  the  explosion 
coupling  into  elastic  waves  - the  feature  presiimably  most 
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indicative  of  explosion  yield.  Evaluation  of  the  comparison 
should  rely  most  heavily  on  the  b phase  data. 

The  comparison  between  theoretical  and  observed  ampli- 
tudes from  the  seismogrcims  of  Figure  3.4  is  shown  in  Figures 
3.5  and  3.6.  The  data  in  Figure  3.5  is  for  the  "b"  phase; 
that  in  Figure  3.6  is  for  the  maximum  aimplitude  in  the  first 
three  cycles,  the  "d"  phase.  Along  with  the  plotted  b and 
d amplitudes,  the  periods  of  the  cycles  at  which  these  were 
measured  appear  on  the  figures. 

With  regard  to  the  data  of  Figures  3.5  and  3.6  there 
are  two  important  observations  to  be  made.  First,  the  ampli- 
tude and  period  of  the  "apparent  ground  motion”  is  dependent 
on  the  characteristics  of  the  recording  instrument.  Unfor- 
tunately, it  is  not  necessarily  true  that  a plot  of  actual 
versus  theoretical  "true  ground  motion"  amplitudes  would 
have  the  same  appearance  as  the  plots  of  Figures  3.5  and  3.6. 
This  is  because  the  instrument  correction  to  true  groxmd 
motion,  as  conventionally  used,  has  an  erratic  effect.  This 
subject  is  discussed  in  some  detail  in  Appendix  A where  we 
show  the  results  of  the  following  nximerical  experiment:  (1) 

The  theoretically  computed  ground  motion  for  KASSERI  at  HNME 
is  filtered  by  six  different  short  period  instruments;  one 
being  specified  by  LRSM  nominal  instrument  response  curves 
and  the  others  specified  by  calibration  data  provided  by  the 
Project  Officer  for  the  five  SDCS  stations.  (2)  Synthetic 
seismograms  were  computed.  (3)  Amplitudes  were  measured  in 
the  conventional  way;  since  the  entire  process  is  done  digit- 
ally, the  errors  in  the  measurement  of  amplitude  and  apparent 
period  can  be  reduced  to  nearly  vanishing. 

The  results  of  the  numerical  experiment  outlined  above 
is  somewhat  surprising.  The  apparent  period  is  used  to  deter- 
mine an  instrument  correction  to  what  we  shall  call  "apparent 
ground  motion" . Even  though  the  actual  ground  motion  is 
identical,  this  apparent  ground  motion  varies  over  a large 
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Comparison  of  theoretical  and  observed  b ampli- 
tudes for  MAST  and  KASSERI  at  five  SDCS  stations 
Seismogrcims  were  computed  for  two  values  of  T/Q 
as  indicated.  The  KASSERI  source  function  is 
denoted  148  and  the  upper  mantle  model  is  HVn^E-3. 
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range;  as  much  as  -a  factor  of  two  1 It  turns  out  that  the 
discrepancy  would  be  much  reduced  if  no  instrument  correc- 
tion were  applied. 

The  second  point  to  be  made  is  closely  related  to  the 
first  made  above.  We  find  that  we  must  know  the  instrument 
response  quite  accurately  to  be  able  to  compute  theoretical 
seismograms  that  are  comparable  to  the  actual  observations. 
But  how  well  do  we  k.icw  this  response?  How  carefully  are 
the  instruments  cali orated?  At  present  we  have  no  way  of 
making  a judgment. 

The  "observed”  data  plotted  in  Figures  3.5  and  3.6 
are  based  on  amplitude  and  period  measurements  made  at  S*. 
Once  again,  the  instrument  correction  problem  arises.  We 
believe  that  the  cimplitudes  appearing  in  the  SDCS  Event 
Reports  (No.  18  for  MAST  and  No.  41  for  KASSERI)  were  cor- 
rected for  instrument  response  using  LRSM  nominal  curves. 

We  used  the  specific  instrument  curves  provided  by  the 
Project  Officer.  These  were  based  on  calibrations  carried 
out  just  before  MAST  and  we  can  only  hope  that  they  remained 
valid  for  KASSERI  - both  for  the  observations  and  for  syn- 
thesized seismograms. 

From  the  comparison  of  Figures  3.5  and  3.5  v/e  draw 
the  following  conclusions; 

• The  scaling  cf  observed  amplitudes  between  the 
two  events  is  quite  consistent  except  for  RKCN. 

We  previously  pointed  out  the  marked  disimilarity 
in  waveforms  at  RKON  for  these  two  events. 

• The  match  between  theoretical  and  observed  ampli- 
tudes is  about  the  scune  for  the  two  events . The 
exceptional  stations  are  RKON  and  FNWV.  It 

will  be  very  difficult  to  explain  the  RKON  anomaly 
with  present  techniques  since  the  epicentral 
distance  variation  between  events  is  only  0.1 
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degrees.  Wavefore  similarities  at  WHY2K,  CPSO 
and  HNME  discourage  attempts  to  attribute  the 
anomaly  to  near-source  effects.  Even  in  our 
theoretical  model  the  interference  pattern  at 
FNWV  is  enough  to  make  the  b cimplitude  behave 
erratically  as  we  see  by  comparing  results  for 
the  two  events.  However,  it  is  unlikely  that 
plane  layered  models  of  the  crustal  structures 
beneath  FNWV  would  permit  a very  accurate  dupli- 
cation of  observed  records  at  this  station. 

The  scaling  between  the  two  computed  source  func 
tions  seems  to  be  approximately  correct.  Per- 
haps the  MAST  source  is  relatively  larger  than 
that  for  KASSERI. 

Differences  between  theoretical  and  observed 
amplitudes  are  less  than  a factor  of  two  nearly 
everywhere.  Minor  adjustments  of  the  upper 
mantle  model  could  improve  the  agreement  at  all 
these  stations. 

Changing  T/Q  from  0.95  to  1.05  increases  the 
period  of  the  phase  measured  by  s 0.03  - 0.07 
seconds.  The  increase  in  amplitude  is  15  - 25 
percent. 

The  periods  of  the  b phase  measured  off  the  syn- 
thetic seismograms  generally  agree  with  the  ob- 
servations. If  anything,  they  may  be  a little 
too  short  on  the  average.  For  the  d phase  the 
theoretical  periods  are  clearly  too  long. 


3.1.5  The  Effect  of  Material  Strength  on  Teleseismic  Body 
wave  Recorainas 
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Tedale  3.1?  We  compare  the  four  sources,  which  vary  only  in 
the  specification  of  the  material  strength,  by  computing 
seismogreims  for  station  HNME  with  T/Q  = 1.05.  The  seismo- 
grams are  shown  in  Figure  3.7.  The  b amplitudes  were  mea- 
sured from  these  records  and  are  plotted  versus  the  limiting 
material  strength  (Y  + Y^)  in  Figure  3.8.  Also  shown  is  the 
period  of  the  b measurement.  Since  these  seismograms  were 
computed  with  the  LRSM  normal  instrument  response,  direct 
comparison  with  the  KASSERI  observations  is  not  possible. 
However,  in  Appendix  A,  Table  A. 2,  we  show  that  the  HMNE 
specific  instrument  gives  a b phase  aunplitude  that  is  0.46 
as  large  and  a period  that  is  0.19  seconds  shorter. 

The  b phase  data  of  Figure  3.8  is  consistent  with  our 
discussion  of  the  spectral  characteristics  of  the  four  source 
functions  (Figure  3.2)  given  in  Section  3.1.1.  The  aunplitude 
increases  with  decreasing  strength  and  the  rate  of  increase 
is  more  rapid  at  lower  strength.  The  period  of  the  waveform 
increases  rapidly  with  decreasing  strength  due  to  the  shift 
in  corner  frequency.  For  b amplitude  there  is  essentially 
no  difference  between  the  two  sources  (146  and  147)  that  have 
the  same  ultimate  strength  as  granite,  even  though  the  shape 
of  the  yield  surface  is  different. 

In  the  previous  section  we  saw  that  148  was  reasonably 
successful.  If  in  error,  one  would  think  148  to  be  too  low 
in  amplitude  and  to  have  too  little  energy  (relatively)  in 
the  high  frequencies.  We  see  in  Figure  3.8  that  improvements 
in  these  two  features  require  adjustments  to  the  strength  in 
opposite  directions.  Also,  a strength  of  about  half  that  of 
granite  seems  reasonable.  If  we  are  satisfied  with  the  MAST 
source  which  scales  to  the  MAST  source  in  a manner  consistent 
with  the  observations  is  also  required.  Source  148  seems  to 
fulfill  this  criterion  fairly  well.  For  all  these  reasons 
mentioned  above,  we  selected  the  source  148  as  being  most 
representative  of  KASSERI.  The  data  of  Figure  3.8  are  quite 


'"■VsriW' 


25 


(rtui)  q 


R-2924 


149 


Figure  3.8.  The  b phase  measurements  for  theoretical  seismo- 
grams at  HNME  computed  with  four  source  functions. 
The  source  functions  represent  KASSERI  tuff  with 
four  different  estimates  of  material  strength. 

The  instrument  response  was  the  LRSM  nominal. 
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important  as  an  indication  of  the  effect  of  erroneous  esti- 
mates of  material  strength. 

3.2  PILEDRIVER  GROUND  MOTION  CALCULATIONS 
3.2.1  Introduction 

An  elastic  wave  propagation  code,  based  on  generalized 
ray  theory,  was  recently  developed.  This  computational  code 
capability  enables  us  to  synthesize  the  complete  elastic  wave 
field  propagating  from  an  explosive  source,  including  both 
near-  and  far-field  terms.  An  immediate  application  of  this 
code  was  made  to  the  PILEDRIVER  explosion,  detonated  in 
granite  at  NTS. 

In  fitting  calculated  ground  motions  produced  by  the 
PILEDRIVER  RDP  to  the  observed  free-field  recordings,  it  was 
assumed  (Cherry,  et  al . , 1975)  that  the  observed  data  were 
relatively  unaffected  by  waves  reflected  from  the  free  sur- 
face or  from  interfaces  between  geologic  layers.  The  RDP 
is,  of  course,  one-dimensional  so  that  such  effects  were  not 
included  in  the  calculations.  The  purpose  of  the  calculation 
to  be  described  below  was  to  ascertain  the  validity  of  this 
assumption.  The  results  of  this  calculation  follow  a brief 
description  of  the  wave  propagation  code. 


3.2.2  Description  of  Wave  Propagation  Code 

The  calculation  of  the  RDP  is  made  with  the  1-D  SKIPPER 
code.  This  code  assumes  that  the  highly  nonlinear  process  in 
which  the  explosion  energy  is  coupled  into  shock  waves  ar I 
then  into  elastic  stress  waves  is  a one-dimensional  process. 
Included  is  the  assumption  that  local  heterogeneities  do  not 
significantly  disturb  the  spherical  symmetry  of  the  radiation. 
It  should  be  noted  that  two-dimensional  sources  may  be  treated 
in  an  analogous  way  (e.g.,  Bache  and  Harkrider,  1976;  Cherry, 
et  al. , 1976) . 
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The  use  of  elastic  wave  propagation  techniques  is 
justified  if  the  following  assumptions  are  correct:  For  the 

wavelengths  of  interest,  the  interfaces  between  layers  are 
plane  parallel.  (This  assumption  may  be  relaxed  to  allow 
non-parallel  layers.)  Although  we  consider  reflections  and 
transmissions  where  the  material  response  is  clearly  in  the 
nonlinear  source  region,  we  assume  that  the  interface  be- 
havior is  not  too  different  from  that  for  elastic  waves  so 
that  elastic  reflection  and  transmission  coefficients  may  be 
used.  The  analytical  propagation  technique  is  exact  insofar 
as  elastic  wave  propagation  in  layered  media  is  concerned. 
Thus,  all  frequencies  are  propagated  accurately  in  the  layered 
medium.  It  is  possible  to  include  an  p.ielastic  attenuation 
factor  in  the  analytical  propagation  theory  (this  is  pre- 
sently being  implemented) , but  no  such  factor  was  included 
in  the  calculations  discussed  here.  Of  course,  the  explo- 
sion wave  is  strongly  attenuated  in  the  1-D  nonlinear  cal- 
culation from  which  the  RDP  is  obtained. 

The  elastic  wave  propagation  program  used  here  is 
based  on  generalized  ray  theory.  A brief  description  of  its 
salient  features  follows.  A ray  theory  is  based  on  the  fact 
that  the  response  of  an  elastic  system  can  be  represented 
as  a series,  each  term  of  which  being  identified  as  the 
contribution  of  a particular  ray.  A ray  is  characterized 
by  a ray  path,  which  is  a description  of  the  nvunber  of  tra- 
verses the  ray  makes  across  the  layers,  and  the  mode  (P  or 

S wave)  of  the  traverses.  For  example,  for  a source  and  re- 

ceiver above  an  interface,  the  solution  is  composed  of  three 
rays:  The  direct  ray,  the  ray  which  travels  to  the  interface 

as  a P wave  and  returns  as  a P wave,  and  the  ray  which  re- 
turns from  the  interface  as  an  S wave.  The  solution  is  formu- 
lated in  a way  described  by  several  authors  (e.g.,  Helmberger, 

1968)  . 

It  has  been  shown  that  the  response  u (a  component 
of  displacement,  velocity,  or  stress)  at  the  cylindrical 
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coordinates  (x,z,9)  and  time  t,  due  to  a source  at  the 
origin  can  be  written  as: 


N 


u(x,z,0;t)  =i/'(t)*^  J f (Pj,T)S(Pj)R(Pj)Tj  (Pj)C(Pj)dT 


t,. 


Here,  j is  the  ray  index,  and  N indicates  the  number  of 
rays  to  be  included.  The  integration  starts  at  the  arrival 
time  tjij*  The  complex  parameter  p^  is  chosen  to  make  the 
imaginary  part  of  t vanish,  where 


i=i  \ Ji  / 


1/2 


(3.4) 


and  and  are  the  layer  velocities  and  thicknesses. 

At  the  arrival  time,  p.(t  .)  is  known  as  the  ray  parameter 

3 0 3 

in  geometrical  optics.  The  function  S(Pj)  specifies  the 
source  type  (e.g.,  explosive),  R(Pj)  specifies  whether  the 
incoming  wave  at  the  receiver  is  shear  or  compressional,  and 
C(Pj)  determines  the  quantity  (e.g.,  vertical  displacement) 
being  computed.  The  effect  of  the  interfaces  is  included  in 
the  product  of  the  reflection  and  tra-.ismission  coefficients 
Tj(Pj),  while  f(Pj,T)  is  a kernal  common  to  all  these  cal- 
culations. The  RDP  (t)  , is  the  time  history  of  the  source 
and  * denotes  convolution. 

The  integration  (3.3)  is  performed  for  each  ray  along 

the  complex  path  set  by  (3.4),  with  the  integrands  computed 

on  a finer  spacing  near  t where  they  are  changing  fastest. 

0 

The  ray  responses  are  then  added  together  and  convolved  with 
ip  to  yield  the  final  result.  Classical  (far-field)  seism- 
ology makes  approximations  to  (3.3),  so  that  the  integration 
need  not  be  performed.  This  is  the  technique  by  which  most 
"ray  theory"  calculations  have  been  done.  These  approximations 
generally  require  that  the  periods  of  interest  be  much  smaller 
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than  the  travel  time,  or  equivalently,  that  the  receiver  be 
many  wavelengths  from  the  source.  For  the  calculations  done 
for  this  report,  the  situation  is  reversed,  making  it  neces- 
sary to  fully  evaluate  (3.3). 

In  evaluating  Eq.  (3.3),  we  must  add  N rays  together 
to  form  the  response.  For  most  layered  geometries,  N is 
infinite.  For  example,  when  a source  and  receiver  are  in  a 
sandwiched  layer,  there  are  an  infinite  number  of  ray  paths 
corresponding  to  all  the  internal  reflections  in  the  layer. 
However,  with  each  reflection  the  rays  are  attenuated  be- 
cause they  travel  farther  and  because  they  leak  energy  out 
of  the  layer.  In  practice,  the  rays  which  are  important  for 
a particular  problem  are  easily  determined.  The  calculation 
is  inexpensive  to  run  so  that  if  one  is  in  doubt  about  the 
inclusion  of  a set  of  rays,  numerical  experiments  can  be 
cheaply  done  to  resolve  that  doubt. 

Note  that  the  RDP  affects  the  solution  via  a convolu- 
tion. Thus,  the  RDP  and  the  layered  responses  are  computed 
independently,  and  finally  convolved. 

3.2.3  Results  of  Source  Calculation 

To  check  whether  reflected  waves  from  the  vertical 
layering,  in  particular  the  free  surface,  could  have  in- 
fluenced free-field  ground  motion  recordings,  the  following 
experiment  was  performed.  Using  the  RDP  as  a source  for  the 
elastic  wave  propagation  code,  the  resulting  wave  field  was 
propagated  in  a layered  medium  to  receivers  at  appropriate 
distance  ranges  from  the  source. 

The  source-receiver  geometry  and  the  elastic  structure 
in  which  they  are  embedded  are  shown  in  Figure  3.9.  The  in- 
formation for  this  figure  is  from  a report  by  Perret  (1968) . 
The  receivers  are  at  stations  B-SL  (range  R = 201  m)  and  16- 
SL  (R  = 470  m) . Comparisons  of  the  radial  motions  resulting 
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(a,S,P)  = (1.439,  .831,  1.52) 


(4.663,  2.71,  2.45) 


Free  Surface 
41  m 

113  m 


(5.852,  3.161,  2.70) 
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R = 201  m 
or  470  m 


460  m 


Figure  3.9,  Piledriver  elastic  layering  sonic-receiver 
geometry.  Numbers  in  parentheses  are  P and 
S velocities  (km/sec)  and  densities  (gr./cmM 
respectively. 
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from  calculations  made  using  the  structure  in  Figure  3.9  with 
those  using  a homogeneous  meditim  are  shown  in  Figure  3.10. 

Note  that  the  free  surface  has  a negligible  effect  on  the 
ground  motion  at  B-SL,  and  a slightly  greater  effect  at  16-SL 
where  the  reflected  ray  is  traveling  more  horizontally  than 
at  B-SL.  Also,  the  RVP's  deduced  from  the  radial  velocities 
in  the  above  comparisons  are  shown  to  be  nearly  the  same 
^Figure  3.11).  Thus,  the  assumption  that  the  ground  motion 
recordings  at  these  stations  are  unaffected  by  geologic  layer- 
ing seems  warranted,  and  the  waveforms  observed  for  PILEDRIVER 
were  appropriate  choices  for  normalizing  the  source  calcula- 
tions as  described  in  Bache,  et  al.  (1975b) . 

3.3  GASBUGGY  SOURCE  CALCULATION 


3.3.1  Introduction 

GASBUGGY,  a nuclear  experiment  designed  to  stimulate 
production  of  natural  gas  from  the  impervious  Pictured  Cliffs 
sandstone  formation,  was  detonated  in  Rio  Arriba  County,  New 
Mexico,  in  a layer  of  Lewis  Shale  (depth  of  burial  of  4240 
feet)  approximately  40  feet  below  the  Pictured  Cliffs  sand- 
stone. The  device  yield  was  29  kilotons.  As  part  of  this 
experiment  a prograim  of  free-field  ground  motion  measure- 
ments was  carried  out  to  define  a seismic  source  function  for 
GASBUGGY.  The  set  of  ground  motion  measurements  obtained 
for  GASBUGGY  (Perret,  1972)  provides  an  excellent  oppor- 
tunity to  further  verify  the  capability  of  the  constitutive 
modeling  embodied  in  the  SKIPPER  code  to  predict  close-in, 
free-field  ground  motion  from  underground  explosions.  The 
results  of  these  types  of  calculations  place  important 
constrains  on  the  source  spectra  that  are  used  for  teleseis- 
mic  ground  motion  predictions. 


Figure  3.11.  Amplitude  spectra  at  two  ranges  comparing  homogeneous  and  layered  cases 
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3.3.2  Source  Calculation 

Free-field  instrumentation,  consisting  of  three- 
component  accelerometers  and  particle  velocity  gauges,  was 
installed  at  four  subsurface  stations,  located  at  slant 
ranges  between  approximately  1500  and  1900  feet  from  the 
GASBUGGY  working  point  (Ferret,  1972).  Figure  3.12  shows 
both  a plan  view  of  the  experimental  arrangement  (left) 
and  a cross-sectional  perspective  of  the  subsurface  station 
locations  and  local  geologic  layering. 

Based  on  the  free-field  ground  motion  measurements. 
Ferret  (1972)  computed  the  reduced  displacement  potential 
(RDF)  for  the  GASBUGGY  source  at  the  four  subsurface  stations 
The  measured  RDFs  are  shown  in  Figure  3.13.  Although  these 
four  stations  were  located  in  different  rock  layers  (Figure 
3.11),  the  residual  values  of  the  RDFs  are  quite  similar; 
differing  from  an  average  value  of  6470  m^  by  at  most  12.5 
percent.  Feak  values  of  RDF  at  these  stations  also  are  rela- 
tively uniform.  This  indicates  that  the  two  dimensional 
effects  of  the  layering  on  RDF  are  quite  small.  Therefore, 
the  GASBUGGY  environment  was  modeled  by  a spherically  sym- 
metric source  calculation  in  an  infinite  layer  of  Lewis 
Shale . 

Some  material  properties  data  for  both  Lewis  Shale 
and  Fictured  Cliffs  sandstone  are  available  from  Terhune 
and  Shaw  (1972).  Figure  3.14  shows  the  failure  envelopes 
for  the  materials  of  interest  and  Figure  3.15  the  compres- 
sibility curves.  A Tillotson  equation  of  state  was  fitted 
to  the  compressibility  curve  for  Lewis  Shale.  Other  rele- 
vant material  properties  are  given  in  Table  3.5.  The  consti- 
tutive modeling,  tensile  failure,  effective  stress  law,  etc. 
were  the  same  as  discussed  for  the  KASSERI  event. 

Figure  3.16  shows  the  RDF  obtained  from  a SKIFFER  cal- 
culation using  the  above  material  properties.  The  RDF  is  in 
excellent  agreement  with  the  measured  values  of  Figure  3.13. 
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TABLE  3.5 

MATERIAL  PROPERTIES  FOR  LEWIS  SHALE 


Bulk  Density 
Air-Filled  Voids 
Sonic  Velocity 
Overburden  Pressure 
Bulk  Modulus 
Shear  Modulus 

Failure  Surface  / Y 

I 'I 


2.48  gms/cm® 
0 

4 . 0 km/sec 
300  bars 
200  kbars 
147.6  kbars 
0 . 2 kbars 

2.0  kbars 

7.0  kbars 
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Subsurface  Station  Locations  and  Geolo; 
instrument  station  locations. 


Reduced  displacement  potential 


Figure  3.13.  Reduced  displacement  potential  records 
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Figure  3.14.  Failure  envelopes  for  GASBUGGY  and  RULISON 
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Figure  3.15.  Compressibility  curves  for  RULISON  and  GASBUGGY 
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Figure  3.17  shows  the  RDP  obtained  using  Pictured  Cliffs 
sandstone  material  properties.  The  high  strength  of  this 
rock  results  in  an  RDP  much  lower  than  measured. 


i 


\ • 


R-2924 


t 


IV.  CONCLUSIONS 

Considerable  progress  was  made  during  the  fourth 
three-month  period  of  this  project  on  several  of  the  tasks 
specified  in  the  work  statement  of  this  contract. 

The  theoretically  computed  and  observed  cimplitudes 
of  both  the  b and  d (maximum)  body  wave  phases  for  KASSERI 
agree  to  well  within  a factor  of  two  at  all  of  the  five 
sues  stations.  Minor  adjustments  of  the  upper  mantle  model 
could  improve  the  agreement  of  all  the  SDCS  sites. 

Computations  of  the  effect  of  material  strength  on 
teleseismic  body  wave  amplitudes  indicate  that  the  amplitude 
of  the  b phase  increases  with  decreasing  strength;  the  rate 
of  increase  being  more  rapid  at  lower  levels  of  material 
strength.  An  additional  effect  is  that  the  apparent  period 
of  the  b phase  increases  rapidly  with  decreasing  strength 
due  to  a shift  in  the  corner  frequency  of  the  source  spectrum. 

Further  verification  of  the  explosion  source  modeling 
code,  SKIPPER,  was  achieved  by  comparison  of  free  field 
ground  motion  calculations  and  observations  for  the  PILEDRIVER 
and  GASBUGGY  explosions.  In  the  case  of  PILEDRIVER,  the 
exercise  of  a recently  developed  elastic  wave  propagation 
code  verified  that  the  observed  PILEDRIVER  RDP  was  not  af- 
fected by  reflected  waves  from  the  free  surface  or  from  layer 
interfaces  in  the  source  region. 

Excellent  free  field  ground  motion  measurements  and 
material  properties  data  are  available  for  the  GASBUGGY  ex- 
plosion. The  computed  RDP,  employing  these  data,  egrees  quite 
well  with  the  measured  RDP,  increasing  our  confidence  in  the 
constitutive'  modeling  of  the  explosion  source  and  ultimately 
the  teleseismic  ground  motion  predictions  based  on  these 
source  calculations. 
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APPENDIX  A 

EFFECT  OF  INSTRUMENT  RESPONSE  ON  MEASURED  AMPLITUDES 
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Standard  procedure  for  determining  ground  notion  ampli- 
tudes from  seismograph  recordings  requires  correction  for  the 
instrument  response  at  the  apparent  period  of  the  cycle  being 
measured.  There  are  several  possible  sources  of  error  in  this 
procedure ; 

1.  The  signal  is  not  monochromatic  and  a single 
period  correction  factor  can  never  entirely  re- 
move the  instrument  effect. 

2.  The  amplitude  and  phase  response  of  the  instru- 
ment at  the  time  the  signal  is  recorded  may  not 
be  accurately  known. 

3.  It  is  often  difficult  to  accurately  measure  the 
dominant  period  of  the  cycle  being  measured. 

These  are  all  sources  of  error  that  refer  to  a single  record- 
ing by  a single  instrument.  While  everyone  knows  they  exist, 
a quantitative  estimate  of  their  effect  is  not  so  widely  known. 

If  we  are  comparing  measured  amplitudes  from  a series 
of  events  that  have  similar  waveforms  and  are  recorded  by  the 
same  (nominal)  instrument,  then  the  errors  1-3  listed  above 
should  be  normally  distributed  and  can  be  accounted  for  by 
statistical  methods.  Let  us  here  pose  a different  question. 
Assume  that  the  same  ground  motion  is  recorded  by  a number 
of  different  instruments  and  that  the  measurements  are  accu- 
rately made.  After  correcting  for  tlie  instrument  response 
at  the  apparent  period  of  the  cycle  being  measured,  how 
closely  do  the  estimates  of  ground  motion  agree?  In  this 
appendix  we  give  the  results  of  several  numerical  experiments 
that  indicate  that  the  differences  can  be  surprisingly  large. 
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First,  let  us  consider  the  effect  on  body  wave  measure- 
ments. The  instruments  to  be  discussed  are  as  follows.  For 
the  standard  of  comparison  let  us  take  the  LRSM  Benioff  short 
period  instrument  for  which  nominal  response  curves  were 
given  in  all  the  old  SDAC  shot  reports.  We  have  been  inter- 
ested in  computing  the  ground  motion  at  five  SDAC  stations, 
RKON,  CPSO,  WHY2K,  FNWV,  HNME.  Calibration  data  from  June  12, 
1975  was  supplied  by  the  project  officer  for  each  of  these 
stations.  This  was  in  the  form  of  s 5 amplitude  response 
values  and  a series  of  phase  response  values.  The  data  were 
not  entirely  consistent  and  had  to  be  smoothed.  Certainly, 
it  is  valid  to  ask  how  trustworthy  these  calibration  data 
are,  but  it  is  the  best  we  have. 

Let  us  view  the  same  explosion-like  ground  motion 
through  the  six  instruments  which  we  denote  as  LRSM  and  by 
the  names  of  the  five  SDCS  stations.  The  ground  motion  is 
that  computed  for  KASSERI  at  HNME  with  T/Q  = 1.05.  The  re- 
sulting seismograms  are  showj.  in  Figure  A.l.  In  terms  of 
signal  shape,  the  main  difference  is  that  distinguishing  the 
LRSM  record  from  the  other  five.  We  also  notice  that  the 
1 Hz  peak  amplitude  values  given  with  each  record  vary  over 
a substantial  range  (1.02  - 1.36  u)  with  the  LRSM  record 
agaiii  Lt.-ing  exceptional. 

What  are  the  amplitudes  associated  with  the  seismo- 
graims  of  Figure  A.l?  With  these  theoretical  records  we  can 
reduce  the  measurement  errors  to  nearly  vanishing.  Rather 
than  measuring  by  eye,  a parabola  is  fit  to  the  digital  data 
defining  the  peaks.  The  amplitude  and  apparent  period  of 
the  cycle  of  interest  can  then  be  determined  with  sources  of 
human  error  removed.  Since  the  instrument  is  input  to  the 
code  in  digital  form,  we  can  determine  the  (theoretical) 
instrument  response  at  the  apparent  period  with  essentially 
zero  error.  The  amplitude  and  period  data  for  the  b and 
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Figure  A.l.  Theoretical  ground  motion  for 
viewed  with  six  different  ins 
strument  designations  appear 
as  I/name.  For  the  calculatio 
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d phases  is  summarized  in  Table  A.l.  The  maximum  or  d phase 
is  indicated  on  the  seismogrcuns  by  a bar. 


Table  A.l.  Amplitude  and  Period  Data  for  the  Records  of 
Figure  A.l 


Instrument 

!b 

b (my) 

d (my) 

I/LRSM 

1.16 

1079 

1.52 

3373 

I/RKON 

0.97 

493 

1.35 

2589 

I /CP SO 

1.05 

607 

1.43 

3093 

I/WiY2K 

0.99 

537 

1.37 

2670 

I/FNWV 

1.06 

686 

1.39 

2725 

I/HNME 

0.97 

491 

l."2 

2390 

The  amplitude  and  period  differences  are  startling. 

If  we  just  consider  the  SDCS  instruments,  the  period  varies 
by  0.1  second  and  the  amplitude  by  ~ 40  percent.  However, 
the  LRSM  gives  an  amplitude  that  can  differ  by  more  than  a 
factor  of  two!  Unfortunately,  the  instrument  correction 
acts  in  the  wrong  way.  If  we  were  to  make  no  instrument 
correction  other  than  to  the  gain  at  1 Hz,  the  discrepancy 
in  amplitudes  in  the  data  of  Table  A.l  would  be  much  reduced. 
For  example,  with  no  correction  the  b amplitudes  for  the 
SDCS  instruments  vary  from  523  - 695  my  and  the  LRSM  ampli- 
tude is  703  my.  Perhaps  we  should  not  make  instrument  cor- 
rections ! 

What  about  when  the  actual  ground  motion  is  different 
from  that  used  in  the  comparison  above?  Let  us  consider 
KASSERI  as  computed  at  each  of  the  SDCS  stations  with  earth 
model  HWNE-3  and  T/Q  = 1.05.  We  compute  seismograms  at  each 
of  the  stations,  first  with  the  LRSM  nominal  instrument, 
then  the  instrument  response  specific  to  each  station.  The 
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results  are  shown  in  Figure  A. 2 and  the  amplitudes  are  sum- 
marized in  Table  A. 2. 


Table  A. 2 

. Amplitude  . 

and  Period 

Data  for 

the  j 

Records  oJ 

Station 

Figure  A. 2 
Instrument 

^b 

b 

^d 

d 

- 

■ 

RKON 

LRSM  Nominal 

1.09 

2682 

1.48 

9345 

Specific 

0.95 

1337 

1.35 

2589 

CPSO 

LRSM  Nominal 

1.29 

2061 

1.84 

6986 

Specific 

1.17 

1249 

1.43 

3093 

WHY2K 

LRSM  Nominal 

1.19 

646 

1.38 

4452 

Specific 

1.03 

336 

1.37 

2670 

FNWV 

LRSM  Nominal 

1.43 

1078 

1.49 

2024 

Specific 

1.38 

732 

1.34 

2725 

HNME 

LRSM  Nominal 

1.16 

1079 

1.52 

3373 

Specific 

0.97 

491 

1.32 

2390 

In 

Table  A. 3 we 

show  the 

same  kind 

of  comparison 

ground  motion  from  the  MAST  event.  In  this  case  the  dominant 
period  of  the  actual  ground  motion  is  somewhat  shorter  than 
for  KASSERI.  Only  the  data  for  the  b amplitude  is  shown  as 
it  illustrates  our  point  — the  instrument  response  can  make 
a big  difference  in  the  apparent  ground  motion  determined 
from  seismograms  by  conventional  means. 


TirE  (SEC) 


I 


R-2924 


TH 


» 


t 


9 


. i 


Table  A. 3.  Amplitude  and  Period  Data  for  MAST  Seismograms 
Computed  with  the  LRSM  Nominal  and  Station 
Specific  Instruments 


Station 

Instrument 

!b 

b 

RKON 

LRSM  Nominal 

0.93 

1154 

Specific 

0.86 

775 

CPSO 

LRSM  Nominal 

0.94 

603 

Specific 

0.92 

478 

WHY2K 

LRSM  Nominal 

0.94 

237 

Specific 

0.89 

167 

FNWV 

LRSM  Nominal 

1.08 

248 

Specific 

1.02 

168 

HNME 

LRSM  Nominal 

0.93 

415 

Specific 

0.86 

270 
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